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Abstract 

This paper summarizes the Multiple Parton Interactions studies in 
CMS, focusing on the already performed low pT QCD measurements 
up to highest centre of mass energies of 7 TeV and discussing the plans 
for the direct measurement of the multiple high-pr scatterings. 

The underlying event in pp interactions is studied measuring the 
charged multiplicity density and the charged energy density in the 
transverse region, which is defined considering the azimuthal distance 
of the reconstructed tracks with respect to the leading track-jet of the 
event, defined from tracks according to a jet clustering algorithms. 

In addition, we present the measurement of the underlying event 
using the jet- area/median approach, demonstrating its sensitivity to 
different underlying event scenarios. 

Observations in the central region are complemented by the mea- 
surement of the energy flow in the forward direction for minimum bias 
and central di-jet events. 

We compare our underlying event and forward results with the 
predictions from different Monte Carlo event generators and tunes, 
aiming to best parametrize the multiple parton interaction energy de- 
pendence starting from the Monte Carlo tunes developed to best fit 
the charged particle spectra measured at central rapidities. 

Finally we discuss the strategy to directly measure the multiple 
particle interactions rate focusing on the topologies with two hard 
scatterings in the same event. 
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1 Introduction 



The Multiple Parton Interactions (MPI) are currently invoked to account 
for observations at hadron colliders that would not be explained otherwise: 
the activity of the Underlying Event (UE), the cross sections of multi-jet 
production, the survival probability of large rapidity gaps in hard diffraction, 
etc. [1]. At the same time, the implementation of the MPI effects in the 
Monte Carlo models is quickly proceeding through an increasing level of 
sophistication and complexity that achieves deep general implications for 
the LHC physics. 

This contribution summarizes the early Underlying Event (UE) and for- 
ward measurements of the CMS collaboration using pp collision data up to 
highest energies of y/s = 7 TeV. It also reports along the feasibility study for 
the direct measurement of double parton scattering phenomena focusing on 
the 3jet + 7 channel. 

A detailed description of the CMS detector is available in Ref. p]. Gen- 
erator level Monte Carlo (MC) predictions are compared to the data cor- 
rected with a bayesian unfolding technique taking into account the detector 
effects [3]. 

The predictions for inelastic events are provided here by several tunes of 
the PYTHIA program, versions 6.420 jH and 8.13^ [7|. PHOJET [H] 
is also used in the forward measurements: 

The pre-LHC tune D6T [9l |10] of PYTHIA-6, which describes the lower 
energy UA5 and Tevatron data, is a widely used reference that will also be 
used for most of the presented analyses. The tunes DW [10] and CW jlT] . 
which were found to describe best the UE CMS data at 0.9 TeV whereas D6T 
predictions were too low [H] , will also be discussed for the 7 TeV data. The 
pre-LHC tune Perugia-0 [12] and the new PYTHIA-6 tune, Zl [13], includes 
Pt ordering of parton showers and the new PYTHIA MPI model [H]. It 
implements the results of the Professor tunes [15j considering the fragmen- 
tation and the color reconnection parameters of the AMBTl tune [TB], while 
the CMS UE results [TT], [T7] discussed in this paper have been used to tune 
the parameters governing the value and the ^/s dependence of the transverse 
momentum cut-off that in PYTHIA regularizes the divergence of the leading 
order scattering amplitude as the final state parton transverse momentum 
Pt approaches 0. The tune Z2 is similar to Zl, except for the transverse 
momentum cut-off at the nominal energy of = 1.8 TeV which is de- 
creased of 0.1 GeV/c. PYTHIA-8 also uses the new PYTHIA MPI model, 
which is interleaved with parton showering. The default Tune 1 (indicated 

ipYXHIA version 8.108 is used in the feasibility studies reported in sectionEl 
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as PYTHIA-8) in and the new PYTHIA-8 tune 4C [E], which focuses on 
the description of the early LHC data, are adopted here. It includes soft and 
hard diffraction (TH], whereas only soft diffraction is included in PYTHIA-6; 
the diffractive contributions are, however, heavily suppressed by the trigger 
and event selection requirements, especially for large pt values of the leading 
track-jet. The PYTHlA-8 tune 4C also focuses on the description of the 
early LHC data. The parton distribution functions used to describe the in- 
teracting protons are the CTEQ6LL set for D6T, Z2 and 4C. The CTEQ5L 
set is adopted for the other simulations [201 1^ . 

2 The Early Underlying Event Measurements 

In the presence of a hard process, characterized by particles or clusters of 
particles with a large transverse momentum px with respect to the beam di- 
rection, the final state of hadron-hadron interactions can be described as the 
superposition of several contributions: products of the partonic hard scatter- 
ing with the highest pt, including initial and final state radiation; hadrons 
produced in additional MPI; "beam-beam remnants" (BBR) resulting from 
the hadronization of the partonic constituents that did not participate in 
other scatterings. Products of MPI and BBR form the UE, which cannot be 
separated from initial and final state radiation. 

The early CMS UE measurements focus on the understanding of the 
UE dynamics studying charged particle production with two different ap- 
proaches. The first (traditional) approach pTl [T7] concentrates on the study 
of the transverse region, which is defined considering the azimuthal distance 
of the reconstructed tracks with respect to the leading track or leading track- 
jet of the event: 60° < |A0| < 120°. The jet reconstruction algorithm used 
in these studies is SisCone [22]. On top of the traditional approach, a new 
methodology using anti- i^T jets [23] and relying on the measurement of their 
area [21] is adopted for the first time by CMS using charged particles in pp 
collision data collected at -y/s = 0.9 TeV The new set of UE observables 
consider the whole pseudorapidity-azimuth plane instead of the transverse 
region and inherently take into account the leading jets of an event. 

The centre-of-mass energy dependence of the hadronic activity in the 
transverse region is presented on the two top Figures [T] as a function of the 
Pt of the leading track-jet. The data points represent the average multiplicity 
and average scalar track-p^ sum dependence, for y/s = 0.9 TeV and a/s = 7 
TeV using tracks with a pseudorapidity \ri\ < 2.0 and px > 0.5 GeV/c. 
A significant growth of the average multiplicity and of the average scalar 
Pt sum of charged particles transverse to that of the leading track-jet is 
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Figure 1: (Upper plots) average multiplicity and average scalar J2Pt in 
the transverse region as a function of the leading track-jet pt, for data at 
^/s = 0.9 TeV and ^/s = 7 TeV. (Bottom left plot) normalized scalar J2pt 
distribution in the transverse region for data at ^/s = 0.9 TeV and ^/s = 7 
TeV; the leading track-jet is required to have pt > S GeV/c. Predictions 
from PYTHIA-6 tune Zl and PYTHIA-8. 135 Tune 4C are compared to the 
corrected data. The inner error bars indicate the statistical uncertainties 
affecting the measurements, the outer error bars thus represent the statistical 
uncertainties on the measurements and the systematic uncertainty affecting 
the MC predictions added in quadrature. (Bottom right plot) normalized 
median of pt over jet area for track-jets reconstructed from collision data at 
^/s = 0.9 TeV (black circles). Predictions from several PYTHIA-6 tunes and 
PYTHIA-8 Tune 1 are compared to data. 
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observed with increasing scale provided by the leading track-jet pt, followed 
by saturation at large values of the scale (more evident for multiplicity profile 
than average scalar sum). A significant growth of the activity in the 
transverse region is also observed, for the same value of the leading track-jet 
Pt, from a/s = 0.9 TeV to a/s = 7 TeV. These observations are consistent 
with the ones obtained at Tevatron [26j. The evolution with the hard scale 
of the ratio of the UE activity at 7 TeV and 0.9 TeV is remarkably well 
described by the Zl MC. The trend is also reproduced by PYTHIA-8 Tune 
4C. The Z2 predictions at y/s = 0.9 TeV (not shown here) agree with Zl in 
shape, but the normalization is 5-10% higher for both the observables; this 
trend is opposite with respect to the one observed at 7 TeV and indicates 
that a less pronounced ^/s dependence of the transverse momentum cut-off 
should be adopted for tunes using the CTEQ6LL PDF set than for the tunes 
optimized for the CTEQ5L set. 

The strong growth of UE activity with charged particles is also striking in 
the comparison of the normalized distributions of charged particle multiplic- 
ity (not shown here) and of scalar pt sum which is presented in bottom-left 
plot of Figure [l] for events at ^/s = 0.9 TeV and ^/s = 7 TeV with leading 
track-jet pt > S GeV/c. The particle pt spectrum (not shown) extends up 
to Pt > 10 GeV/c, indicating the presence of a hard component in particle 
production in the transverse region. The distributions for track-jet > 3 
GeV/c, which extend up to quite large values of the selected observables 
in the transverse region are quite well described by the various MC mod- 
els, over several orders of magnitude. This observation gives support to the 
implementation of MPI in PYTHIA. 

The novel technique to quote the UE activity [27] relies on the introduc- 
tion of "ghosts" , virtual deposits of very low energy filling the overall phase 
space which are taken into account by the jet clustering algorithm. The es- 
timator of the overall soft background activity in an event can be derived as 
the median of the ratio between the transverse momentum and the area of 
the jets. One of the advantages of using the median compared to the mean 
is that it turns out to be less sensitive to the influence of outliers, i.e. in par- 
ticular the leading jets in an event. In order to cope with the low occupancy 
observed at y/s = 0.9 TeV, CMS redefines such observable restricting the 
median only to those jets which have physical deposits on top of ghosts: 



p' = median 



C (1) 



A, 

where C is the occupancy of the event, which is the summed area 
covered these jets divided by the considered detector region Atot- In the 
CMS analysis at = 0.9 TeV, jets are reconstructed with the anti- i^T 
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algorithm using tracks with Ir^l < 2.0 and pt > 0.3 GeV/c. In the right 
bottom plot of Figure [T] the p' observable is presented for minimum bias 
events. The general pattern of deviations from data with respect to the 
considered PYTHIA tunes looks rather similar to the one observed with the 
traditional UE measurement. 

3 Study of the Activity in the Forward Re- 
gion 

CMS reports a measurement of the energy flow in the forward region (3.15 < 
|?7| < 4.9, where t] denotes the pseudorapidity) [28] for minimum bias and 
dijet events in pp interactions with centre-of-mass energies ^/s of 0.9 TeV, 
2.36 TeV and 7 TeV. This measurement is connected to the ones reported in 
the previous sections as the basic philosophy is the same: it concentrates on 
the complementary activity of a pp interaction for different energy scales of 
the reconstructed leading objects. 

The energy flow in the region of the Hadron Forward detector is measured 
in two different event classes: in minimum bias events and in events with a 
hard scale provided by a dijet system at central pseudorapidities (|?7| < 2.5) 
and with transverse energy Exjet > 8 GeV for y/s = 0.9 TeV and 2.36 TeV; 
the dijet threshold is increased to 20 GeV for ^/s = 7 TeV. The results are 
qualitatively similar at all the investigated centre of mass energies. Fig. [2] 
shows the results of the forward energy flow at ^/s = 7 TeV for the two event 
classes compared to predictions from Monte Carlo event generators. The 
measured forward energy flow is found to be significantly different between 
the two event classes, with a sensitive increase and a more central activity 
seen in dijet events. 

4 Multiple Parton Interactions in High-pT Phe- 
nomenology 

Quantifying the MPI cross sections basically deals with the measurement of 
(Jeff, the scale factor which characterizes the inclusive rate of the interac- 
tions [291 EO]- From a phenomenological point of view cTe// is a non per- 
turbative quantity related to the transverse size of the hadrons and has the 
dimensions of a cross section. The measurements performed by the AFS, 
CDF and DO collaborations [311 E21 E31 EH ES] favor smaller values of CTg// 
with respect to the naive expectations. The consequent increased rates of 
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Figure 2: Energy flow in the minimum bias (left) and di-jet (right) samples 
as a function of rjobs at a/s = 7 TeV. Uncorrected data are shown as points, 
the histograms correspond to the MC predictions. Error bars corresponds to 
statistical errors. The shaded bands in these plots represent the systematic 
uncertainties of the measurements, which are largely correlated point-to- 
point. 



multiple parton interactions can be interpreted as an effect of the hadron 
structure in transverse plane [3^. Extending such measurements at the LHC 
and studying the possible scale dependency of cXe// is definitely of great in- 
terest and may have a deep impact on the data driven estimations of the 
MPI backgrounds to new physics. 

The production of four high-p-r jets is the most prominent process to 
search for multiple high pt scatterings: two independent scatters in the same 
pp or pp collision {double-parton scattering, DPS) each producing two jets. 
Such a signature has been searched for by the AFS experiment at the CERN 
ISR, by the UA2 experiment at the CERN SppS and by the CDF and DO 
experiment at the Fermilab Tevatron. 

However, searches for double-parton scattering in four-jet events at hadron 
colliders may face significant backgrounds from other sources of jet produc- 
tion, in particular from QCD bremsstrahlung (Fig. |3]-left). Typical thresh- 
olds employed in jet triggers bias the event sample towards hard scatterings. 
However, a high-p^ jet parton is more likely to radiate additional partons, 
thus producing further jets. Thus, the relative fraction of jets from final-state 
showers above a given threshold is enlarged in jet trigger streams which is an 
unwanted bias. On the other hand, looking for four jets in a minimum-bias 
stream will yield little statistics. 

Therefore the strategy to directly measure the MPI rate in high-p^ regime 
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Figure 3: Definition of azimutlial angle between pairs, togetlier with typical 
configurations of double-bremsstrahlung (left) and double-parton scattering 
events (right). 



at hadron colliders also includes the study of multi-jets or jets+photon fi- 
nal states. Indeed the CDF and DO collaborations studied final states with 
one photon and three jets looking for pairwise balanced photon-jet and dijet 
combinations. The data sample was selected with the experiment's inclusive 
photon trigger, thereby avoiding a bias on the jet energy. The better en- 
ergy resolution of photons compared to jets purifies the identification of Et 
balanced pairs. Tevatron found an excess in pairs that are uncorrelated in 
azimuth with respect to the predictions from models without additional hard 
parton scatters per proton-proton scatter. CDF interpreted this result as an 
observation of double-parton-scatters. 

Analyses trying to identify two hard scatters in multi-jet events typically 
rely on methodologies which overcome combinatorics. There are three pos- 
sible ways to group four objects into two pairs: combinations are commonly 
selected pairing objects which are balanced in azimuth and energy. The flavor 
or other specific features of the jets may be used to decrease the combina- 
torics and to make looser the constraints on the balancing. One example of 
such a final state is constituted by events with two h jets and two additional 
light jets. 

In order to discriminate double-parton scatters against double-brems- 
strahlung events, CMS studies the observables A0(~\ employed by AFS, and 
A0*^+\ employed by CDF, probing the azimuthal angle between pairs (Fig.[3|. 
Expectations for the above described variables are therefore A0*^~) -k /2 and 
A^*^"*") TT if additional jets come from double-bremsstrahlung. Otherwise, 
i. e. if additional jets come from multiple interactions, both variables should 
be distributed uniformly. 
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Figure 4: Differential cross section shape as a function of A(f)^~^ (upper 
plots) and A^*^"*"^ (bottom plots) variables. Predictions from PYTHIA 8.108 
{Default scenario) and HERWIG 2.2.0 (left panel) and from three different 
Pythia settings (right panel) shown. 
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Differential cross section shape predictions for tlie A0 observables in pp 
interactions at 14 TeV are sfiown in Fig. |4| HERWIG 2.2.0 |371 |38] and 
PYTHIA-8.108 with default settings which include multiple interactions and 
showering predict similar shapes (Fig. |4]-left). The discrimination power of 
the selected observables to Multiple Parton Interaction patterns is clearly 
shown in Fig. |4]- right, where events with MPI switched off (Shower scenario) 
are compared to events with parton shower switched off [MI scenario). The 
differences are particularly pronounced when selecting the A(f)^^^ observable. 



5 Conclusions 



A strong growth of the UE activity is observed with increasing leading track- 
jet pt for both ^/s = 7 TeV and ^/s = 0.9 TeV. At 7 TeV this fast rise 
is followed above ~8 GeV/c by a saturation region with nearly constant 
multiplicity and small px increase. The same pattern is observed at 0.9 
TeV, with the saturation region starting at ~4 GeV/c. A strong growth 
of the activity is also observed with increasing centre-of-mass energy. The 
large increase of activity in the transverse region is also observed in the J2 Pt 
distribution, indicating the presence of a hard component in the transverse 
region. Very good post-LHC MC tunes are available for the description of 
the UE in the central region. 

A measurement of the underlying event using the jet-area/median ap- 
proach is also reported, demonstrating its sensitivity to different underlying 
event scenarios. 

Complementary underlying event measurements in the forward region 
are also presented. The energy flow in the forward direction is measured 
for minimum bias and central di-jets events. A more global UE description 
including both the central and the forward regions is certainly one of the 
next MC tuning challenges, with deep impact on the understanding of the 
MPI dynamics. 

The Multiple Parton Interactions measurement strategy in the high-p^ 
regime is also briefly discussed focusing on the 3jet + 7 topology. The very 
good performances of the LHC machine should allow to have soon the inte- 
grated luminosity conditions adequate to perform these measurements over a 
wide range of energy scales, with deep impact on the data driven estimation 
of the MPI backgrounds to searches. 
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